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■ We discuss the radiation problem of total reflection for a time-harmonic gen- 
1^ , eralized Maxwell system in a nonsmooth exterior domain C , N > 3 , with 

I nonsmooth inhomogeneous, anisotropic coefficients converging near infinity with a 
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exponentially at infinity. We prove that the corresponding eigenvalues do not accu- 
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I to zero. Finally we are able to generalize these results easily to the corresponding 
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1 Introduction 

If we choose a time-harmonic ansatz (resp. Fourier transform with respect to time) for 
the classical time dependent Maxwell system in 

-curlH + 9tD = I , curlE + 9tB = 

div D = p , div B = 

we are led to consider the time-harmonic Maxwell system with non zero complex frequency 
u and complex valued data e , /i , J and p 

-cnilH + iujeE = I , cm\ E + iufiH = , (1.1) 

div eE = p , divfiH = . (1.2) 

This ansatz may be justified by the principle of limiting amplitude introduced by Eidus 
in [3]. Here we denote the electric resp. magnetic field by E resp. H , the displacement 
current resp. magnetic induction hj D = eE resp. B = pH and the current resp. charge 
density by / resp. p . The matrix valued functions e and p are assumed to be time 
independent and describe material properties, i.e. the dielectricity and permeability of 
the medium, curl = Vx (rotation) and div = V ■ (divergence) mark the usual differential 
operators from classical vector analysis. By differentiation we get 

div eE = — - div / , div pH = 

u 

from (11. ip . such that we can neglect (for co 0) the equations (11. 2p . To formulate these 
equations as a boundary value problem in a domain i7 C we need a boundary condition 
at dQ. Modeling total reflection of the electric fleld at the boundary, i.e. \ f2 is a 
perfect conductor, we impose the homogeneous boundary condition (assuming suflicient 
smoothness of the boundary for the purpose of these introductory remarks) 

uxE = on dn , (1.3) 

which means that E possesses vanishing tangential components aX dQ . Here u denotes 
the outward unit normal on dQ and x the vector product in . We are interested in 
the case of an exterior domain Q , i.e. a connected open set with compact complement. 
Therefore we have to impose an additional condition like 



^xH + E,^xE-H = o{r'^) 



(1.4) 
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(^(x) := x/\x\ , r{x) := the classical so called outgoing Silver-Miiller radiation con- 
dition, which allows to separate outgoing from incoming waves. Interchanging + and — 
in (11. 4p would yield incoming waves. We call the problem of finding E and H with (II. ip . 
(II. 3p and (II. 4p the radiation problem of total reflection for the time-harmonic Maxwell 
system. 

In 1952 Hermann Weyl [21] suggests a generalization of the system (11.11) and (11.31) on 
Riemannian manifolds Q of arbitrary dimension with the aid of alternating differential 
forms. If is a form of rank q (g-form) and H a (g-|- l)-form and if we denote the exterior 
differential dresp. the codifferential S (acting on q- resp. (g + l)-forms) by 



rot := d 



resp. 



div:=5 = (-l)'^^*d* 



to remind of the electro-magnetic background (*: Hodge star-operator), the generalization 
of our system (II. ip and (11.31) reads 



div if + i ueE = F 

L*E = 



votE + iufxH = G 



(1.5) 
(1.6) 



and we call it the generalized time-harmonic Maxwell system of total reflection. Now F 
(former J) is a g-form, G (former 0) a (g + l)-form, e resp. /i a linear transformation 
on g- resp. (g + l)-forms, l : dQ ^ Q the natural embedding and l* the pull-back of 
L . In the case = 3 and g = 1 , i.e. ii^ is a 1-form and H a 2-form, the generalized 
Maxwell system is equivalent to the classical Maxwell system of a perfect conductor, since 
the operators rot = d and div = S acting on g-forms are nothing else than the classical 
differential operators curl and div if g = 1 resp. div and — curl if g = 2 . Moreover, for 
= 3 and 1- resp. 2-forms E we observe that the boundary condition (ll.6p means in the 
classical language u x E = resp. u ■ E = on the boundary, i.e. vanishing tangential 
resp. normal components of the considered fields. We remark that another classical case 
is discussed by this generalization. If iV = 3 and g = resp. g = 2 , i.e. E resp. H are 
scalar valued, we get the equations of linear acoustics with homogeneous Dirichlet- resp. 
Neumann boundary condition, because rot = dresp. div = 6 turns out to be the classical 
gradient V on 0- resp. 3-forms. Moreover, rot = d resp. div = 5 is the zero-mapping 
on 3- resp. 0-forms. In the case of an exterior domain Q C , which we want to treat 
in this paper, we give a generalization of the radiation condition (II. 4p later. For a short 
notation we introduce the formal matrix operators 



M :-- 



' div' 
rot 



A 



e 
/i 



acting on pairs of g-(g + l)-forms and write our problem (II. 5p . (II. 6p easily as 
{M + iujA){E,H) = {F,G) , i*E = 



;i.7) 



:ii 



(For typographical reasons we write form-pairs as {E, H) , although the matrix calculus 



would expect the notation 



H 
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Time-harmonic exterior boundary value problems concerning the classical Maxwell 
equations, i.e. N = 3 and q = 1 , have been studied by Miiller [12] in domains with 
smooth boundaries and homogeneous, isotropic media, i.e. e = /i = Id , with integral 
equation methods and by Leis [7j (see also [9j) with the aid of the limiting absorption 
principle for media, which are inhomogeneous and anisotropic within a bounded subset 
of Q . The generalized time-harmonic Maxwell system has been treated by Week [26j and 
Picard [I?]. 

In this paper we want to discuss the time-harmonic radiation boundary value problem 
of total reflection for the generalized Maxwell equations (11. 8p in an exterior domain Q 
of for arbitrary dimensions N and ranks q . A main goal of our investigations is to 
treat data {F,G) in weighted L^(n)-spaces and inhomogeneous, anisotropic and irregular 
(L°°(f2)-) coefficients e , /i converging near infinity with a rate r~'^ , r > , towards the 
identity. (r(x) := |x| denotes the Euclidean norm in .) We follow in close lines the 
papers of Week and Witsch [30] and Picard, Week and Witsch [22l part 1], which deal 
with the system of generalized linear elasticity and the classical Maxwell equations. In 
particular we generalize the results obtained in the second paper to arbitrary dimensions 
N and ranks of forms q . To present a time-harmonic solution theory we prove that for 
nonzero frequencies to and data {F,G) E L^'|(n) x L^'|'''^(r2)E and L°°-coefficients e , fx 

^2 ^2 

a Fredholm alternative holds true. The main tool to handle irregular coefficients is a 
decomposition lemma, which allows us to prove the polynomial decay of eigensolutions 
as well as an a-priori estimate needed to establish the validity of the limiting absorption 
principle by reduction to the similar results known for the scalar Helmholtz equation. 
The key to this decomposition lemma are weighted Hodge-Helmholtz decompositions, i.e. 
decompositions in irrotational and solenoidal fields, in the whole space case, which have 
been proved in |29j . 

The idea of the decomposition lemma is to use a well known procedure to decouple 
the electric and magnetic field by discussing a second order elliptic system. To illustrate 
this calculation let us look at (11.81) in the homogeneous case A = Id . Applying M — iu 
yields 

{M^ + uj^){E,H) = (M -iLo){F,G) . (1.9) 

If we choose F solenoidal, i.e. div F = , and G irrotational, i.e. rot G = , these 
properties will be transfered to E , i.e. div E = , and H , i.e. rot H = , by (11.81) 
because of 

div div = and rot rot = 

From A = rot div + div rot , where the Laplacian acts on each Euclidean component, we 
get M'^{E,H) = (div rot rot div ^T) = A{E,H) and finally (KM turns to the (compo- 
nentwise) Helmholtz equation 

{A + uj^){E,H) = {M -iuj){F,G) . (1.10) 

Armed with the polynomial decay of eigensolutions and an a-priori estimate for the solu- 
tions corresponding to non-real frequencies (We get these solutions from the existence of a 



*The Definitions will be supplied in section [5] 
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selfadjoint realization of M .) we obtain our radiating solutions for frequencies u G ]R\{0} 
with the method of limiting absorption invented by Eidus [2] as limits of solutions for fre- 
quencies u G C+ \ M . We have to admit finite dimensional eigenspaces for certain eigen- 
values but show that these possibly existing eigenvalues do not accumulate in M\ {0} . All 
these results can be proved by the techniques used in [22] and for orders of decay r > 1 . 
Thus we do not want to repeat them in this paper. However, we refer the interested 
reader to [131 Kapitel 4] for the detailed proofs. 

Proving an estimate for the solutions of the homogeneous, isotropic whole space prob- 
lem with the aid of a representation formula and studying some special convolution kernels 
(Hankel functions) we even can exclude as an accumulation point of eigenvalues. Thus 
the time-harmonic solution operator is well defined on L^'|(i7) x L^'|^^(fi) for small 

frequencies w 7^ . To reach this aim we have to increase the order of decay of the coef- 
ficients e — Id , /i — Id to r > {N + l)/2 and assume that they are in the outside of 
an arbitrarily large ball. Assuming stronger differentiability assumptions on e — Id and 
fi — ld, i.e. in the outside of a ball, we are able to show the exponential decay of eigen- 
solutions as well. To the best of our knowledge it is an open question whether there exist 
such eigenvalues in this general case. Recently under comparable stronger assumptions 
on the coefficients Bauer p] was able to prove that no eigenvalues occur in the classical 
case of Maxwell equations (A^ = 3 , g = 1). Unfortunately his methods are not applicable 
in our general case. It seems to be the same problem that arises trying to prove the 
principle of unique continuation for the generalized Maxwell equation. In the classical 
case the principle of unique continuation was shown by Leis [S] or [SI p. 168, Theorem 
8.17]. However, in the case of homogeneous, isotropic coefficients, i.e. e = Id, /i = Id, 
in the outside of a ball all components of a possible eigensolution solve the homogeneous 
Helmholtz equation (compare (11.101) ) near infinity and therefore by Rellich's estimate [S] 
must have compact support. With the validity of the principle of unique continuation for 
our Maxwell system this eigensolution must vanish. In the general case the principle of 
unique continuation is valid for scalar valued C^-functions e , fi and in the classical case 
for matrices e , n with entries in . (See the citation above from Leis.) 

Having established the time-harmonic solution theory in section [3] we approach the 
low frequency asymptotics of our time-harmonic solution operator. To this end first we 
have to provide a static solution theory. This one is more complicated than for example 
the static solution theory for Helmholtz' equation. The first reason is that for u = the 
system (11.51) resp. (11.81) . i.e. 

TotE = G , diYH = F 

is no longer coupled and that we have to add two more equations to determine E and H , 
i.e. 

div eE = f , rot iiH = g , (1.11) 

which in the case u ^ automatically follow by differentiation from (11.51) as mentioned 
above. (/ = ~~ div F and 9 = —-^ rot G , if div F and rot G exist.) Furthermore, we need 
a boundary condition for the magnetic field (form). Because rot = dand l* commute we 
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derive iujL*^H = l*G for u ^ from fll.Sp . This suggests to impose a condition on the 
term l*^H and, for example, we can choose the homogeneous boundary condition 

L*fxH = 

for our magnetic field. The second reason is that this static Maxwell boundary value 
problem 

TotE = G , divH = F 

dweE = f , TotfiH = g , (1.12) 

L*E = , L*fxH = 

has a nontrivial kernel ^^'^i^) x Ai^^^-i^^^^(^) consisting of harmonic Dirichlet forms. 
Thus we are forced to work with orthogonality constraints on the static solutions to 
achieve uniqueness. For the static system (11.121) a solution theory was given by Kress [6] 
and Picard [16] for the homogeneous, isotropic case, i.e. e = Id , /i = Id , by Picard [21] 
for the inhomogeneous, anisotropic case (Here e and /i even are allowed to be nonlinear 
transformations.) as well as by Picard [TS] for the inhomogeneous, anisotropic classical 
case. For our purpose we need a result like that given by Picard in [16]. In [Tlj we will 
discuss the electro-magneto static problem with inhomogeneous, anisotropic coefficients 
e , ^ in detail. We shortly present some of these results and introduce our static solution 
concept in section HI 

Then in section El the main section of this paper, we prove the convergence of the 
time-harmonic solutions to a special static solution of (I1.12p . This result generalizes the 
paper of Picard [20j, which considers the classical Maxwell equations, to arbitrary odd 
dimensions and ranks g 7^ as well as to coefficients and right hand side data, which 
necessarily do not have to be compactly supported. We note that similar results hold true 
for even dimensions. Since the complexity of the calculations increases considerably due 
to the appearance of logarithmic terms in the fundamental solution (Hankel's function), 
we restrict our considerations to odd dimensions. 

The last section [6] deals with inhomogeneous boundary conditions. Using a new result 
from Week [28], which allows to define traces of g-forms on domains with Lipschitz- 
boundaries, we discuss the time-harmonic problem 

{M + iujA){E,H) = {F,G) , L*E = X 

and the static problem 

M{E, H) = {F, G) , {i*E, i*ijH) = (A, x) 

idiYeE,TotfiH) = if,g) 

It turns out that the solution theories as well as the low frequency asymptotics for these 
problems are easy consequences of the results for homogeneous boundary conditions and 
the existence of an adequate extension operator for our traces. 
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Easily by the Hodge star-operator we always get the corresponding dual results, but 
we renounce them to shorten this paper. 

Essentially this is the first part of the authors ph. d. thesis. Thus sometimes we only 
sketch or neglect some proofs and do not mention all results obtained in To get more 
details on the proofs or some additional results we refer the interested reader to \13\ . 

This paper is the first one in a series of three papers having the aim to determine 
the low frequency asymptotics of the solutions of the time-harmonic Maxwell equations 
completely. In the second paper [T3] we will discuss the corresponding electro-magneto 
static equations in detail and show, how one may define powers of a static solution op- 
erator in weighted Sobolev spaces. This allows us to write down a generalized Neumann 
sum, which is a good candidate for the asymptotic series approaching the time-harmonic 
solutions for small frequencies. In the third paper we finally present the complete low 
frequency asymptotics in the operator norm of weighted Sobolev spaces up to arbitrary 
orders in powers of the frequency. 



2 Definitions and preliminaries 

We will consider an exterior domain Q C , i.e. \ i7 is compact, as a special 
Riemannian manifold of dimension 3 < G N . 

We fix a radius ro and some radii r„ := 2"'ro , n G N , such that \ i7 is a compact 
subset of , the open ball with radius tq centered at the origin. For later purpose we 
choose a cut-off function rj , such that 

r;GC°°(M,M) , supp r; C [1, oo) , ^l[2,oo) = 1 > (2-1) 
and define two other cut-off functions by 

r/(t):=r;(l + ^^) (2.2) 
r2 — ri ^ 

and 

Tj ■= r) or . (2.3) 

Setting Ar := \ Ur and Zr^f := Ar H Uf we note supp Vrj C Zr-^^r2 ■ 
Using the weight function 

p:=(l + r2)V2 

we introduce for m G No and s G M the weighted Sobolev spaces 

H^(fi) := {u G Ll^{n) : p'+l"' 9" m G L\Q) for all |a| < m} 

C H^(fi) := {u G LL(fi) •.p'd'^ue L\Q) for all |a| < m} 

Equipped with their natural norms these are clearly Hilbert spaces. In the special cases 
m = or s = we also write 

H™(fi) := H™(n) , H™(fi) = H™(fi) 
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In Q we have a global chart, the identity, and thus naturally Q becomes a A^-dimensional 
smooth Riemannian manifold with Cartesian coordinates {xi, . . . ,Xn} . For alternating 
differential forms of rank g G Z (g-forms) we define componentwise partial derivatives 
(9" M = {d" ui) dx^ , if M = dx^ (sum convention!), where / are ordered multi-indices of 
length q , and introduce for m e Nq and s e M the weighted Sobolev spaces H™''^(f2) resp. 
H™''?(f2) of g- forms. (Clearly we use the natural componentwise norms in these Hilbert 
spaces.) Again in the cases m = or s = we use the same abbreviations as in the scalar 
case. Especially for m = s = and / = fidx^,g — gidx^ e L^'5(n) we have the scalar 
product 



= / f^*9= / *{f,9)q= / {f:9)qdX= / flQjdX 

Jn Jci Jn Jn 

(A : Lebesgue-measure, ( • , • )g : pointwise scalar product, * : Hodge star-operator) 

Throughout this paper wc denote the exterior derivative dby rot and the co-derivative 
5 — ± * d* hj div to remind of the electro-magnetic background. Because of Stokes' 

o 

theorem and the product rule on C^'^(f2) (the vector space of all smooth g-forms with 
compact support in flj these linear operators are formally skew adjoint to each other, i.e. 

(rot$,^)L2,.+i(Q) = -(<l>,div^)L2.,(n) V $,^'eC~'«(fi) , (2.4) 

which gives rise to weak definitions of rot and div . We note that still rot rot = , 
div div = and rot div -|- div rot = A hold true in the weak sense. Furthermore, for s G M 
we need some special weighted spaces suited for Maxwell's equations: 



Km 

c Rf(0) 



= {EeLl'^{n):TotEeLltt\^)} 
^ {Ee L2'«(Q) : rot£; G LI'1+\Q)} 

^{HeLl'%n): divHeLlf^\n)} 



Equipped with their natural graph norms these are all Hilbert spaces. To generalize 

o o 

the homogeneous boundary condition wc introduce R^(^^) resp. Rf(^^) as the closure of 

o 

C°°''^(0) in the corresponding graph norm || • resp. || • || . Using Stokes' theorem 

we see that in fact the homogeneous boundary condition = is generalized in these 

o 

spaces. The spaces Rs(^^) , Ds(f^) and even Rs(f2) are invariant under multiplication 
with bounded smooth functions (/?, i.e. for E G Rf(fi) we compute 

rot{(pE) = (rot (p) AE + ipiotE 

A subscript at the lower left corner indicates vanishing rotation resp. divergence, e.g. 

o o 

oR^(il) = {E G R^(S1) : rot£^ = 0} , and in the special case s = we neglect the weight 
index, e.g. oD^(^) oDo(^) ■ '^^ consider the whole space, i.e. Q — , we omit the 
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dependence on the domain and write for example qR-s := oR-Kl^^) • For every weighted 
Sobolev spaces Vt , t eM., we define 

t<s t>s 

e.g. D'' 1 (1^) . Moreover, replacing a weight index s by the symbol s = loc resp. s = vox 

means that the forms are locally square-integrable at infinity but square-integrable up to 
the boundary resp. that the forms have bounded supports. 

Furthermore, the rule of partial integration (12. 4p may be generalized as follows: Using 
a usual cutting technique we get for 

{E, H) e R^tiQ) X resp. {E, H) e R?(f]) x Df ^(r^) 

with t, s G M and t + s > resp. t + s > —1 

(rot E, i/)L2.«+i(n) + {E, div iy)L2,,(n) = . (2.5) 
Now let us introduce our transformations: 

Definition 2.1 Let r > . We call a transformation e r-admissible, if 

• e{x) is a linear transformation on q-forms for all x eVL , 

• e possesses L°°{Q)- coefficients, i.e. the matrix representation of e corresponding to 
the canonical basis {^and then for every chart basis {d/i^}) has L°°{Q)- entries, 

• e is symmetric, i.e. for all E, H E L^'''(fi) 

(erE, if)L2,9(Q) = {E,eH)i,2,q(^Q) 
holds, and uniformly positive definite, i. e. 

V A (^-^'-^)L2.9(n) > C- ||-E||l2,,(q) , 

c>o seL2.'j(n) 

• e is asymptotically the identity, i.e. e = eold+e with Eq G M+ and e = 0(r^^) as 
r — )■ oo . We call t the 'order of decay ' of the perturbation e . 

For some results obtained in this paper we need one more additional assumption on 
the perturbations e of our transformations. That is e has to be differentiable in the outside 
of an arbitrarily large ball. More precisely: 

Definition 2.2 Let r > . We call a transformation e r-C^-admissible, if 

• e is T-admissible 
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• and e G C^{Arg) , which means that the matrix representation of i corresponding to 
the canonical basis (^and then for every chart basis {d/i^}) has { A r^) -entries, with 
the additional asymptotic 

dni = 0{r~^"') as r oo , n = l,...,N 

Moreover, we need a special property of our boundary d Q : 

Definition 2.3 A bounded domain S possesses the 'Maxwell compactness property', 
shortly MCP, if and only if the embeddings 

are compact for all q . 

The MCP is a property of the boundary and there is a great amount of hterature about 
the MCP. The first idea was to estimate the H-'^'''(S)-norm by the (R^(S) fl D'^(S))-norm 
(Gaffney's inequahty) and then to use Relhch's selection theorem. To do this one needs 
smooth boundaries, which, for instance, may be seen in P, p. 157, Theorem 8.6]. If g = 
we even have 

R°(S) n D°(S) = R°(S) = Hi'°(H) 

In 1972 [261 127] Week presented for the first time a proof of the MCP for bounded manifolds 
with nonsmooth boundaries ('cone-property'). More proofs of the MCP were given by 
Picard [T9j ('Lipschitz-do mains') and in the classical case by Weber [25] (another 'cone- 
property') and Witsch ^ ('p-cusp-property). A proof of the MCP in the classical case 
for bounded domains handling the largest known class of boundaries has been given by 
Picard, Week and Witsch in [22]. They combined the techniques from [27 [ fTOt 132]. 

Definition 2.4 VL possesses the 'Maxwell local compactness property', shortly MLCP, if 
and only if the embeddings 

are compact for all q . 

Remark 2.5 The following assertions are equivalent: 
(i) possesses the MLCP. 

(ii) Vtr\Ug possesses the MCP for all g > . 

(iii) The embeddings 

are compact for all t, s G M with t < s and all q . 
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(iv) For all t, s G M with t < s , all q and all 0-admissible Sq the embeddings 
are compact. 

Let e be a 0-admissible transformation and t G R. We introduce the ^(weighted 
harmonic) Dirichlet forms' 

e^K^n) ■.= okmne-\Btin) (2.6) 

and in the special case e = Id we denote them by ^K^(fi) . If t = , we always write 
.^K'^in) := ,:K^(fi) . 

o 

By the projection theorem and the L^'^(r2)-orthogonahty of rotR'?~^(i7) and oD^{n) 

o o o 

resp. divD'?+i(i7) and oR''(fi) as well as the inclusions rotR'?-i(i7) c oR-^(^) and 
divD'?+^(i7) C oD^(fi) we get the following Helmholtz decompositions: 

L2''?(fi) = rot R'^-l(^]) e-\D\Q) = oR''{Q) e-^liv 

= e-hot R'?-l(^]) ®e oD'(fi) = e-\R''{n) div D9+i(l]) 

5 ^'^■^^ 

= rot R9-i(n) e^'i^) ®e e'^div D«+i(f]) 

= £-^rot R'?-l(^]) ©, e-\-i:K'^{n) ®s div D'?+i(r]) 

Here all closures are taken in L^''^(fi) and we denote the {e ■ , ■ )L2,<j(Q)-orthogonality by ©^ 
and put © := (Bid ■ These Helmholtz decompositions may be found in [16^ Lemma 1], pT| 
Lemma 1] or in the classical case in [HI p. 168], |22l Lemma 3.13]. 

If Q possesses the MLCP and e is r-C^-admissible with r > , then [131 Lemma 3.8] 
shows 

,:K%(^]) = = s^InM) (2.8) 

and using the Helmholtz decompositions fl2.7p we easily see 

ci"? := dim,:K''(fi) = dim:K5(fi) < oo 

i.e. depends neither on weights —N/2 < t < N/2 — 1 nor on the transformation e . 
Finally we define three operators 

i?:=rdr A ■ =x„dx" A ■ , T:={-iy^*R* , S : = 

acting pointwise on q- resp. (g + 1)- resp. pairs of q- and {q + l)-forms, which will be 
useful to formulate the radiation condition. These operators correspond to rot , div and 



T 
R 



(2.9) 
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M in the following way: If is a smooth function, E a g-form with weak rotation and H 
a (g + l)-form with weak divergence, then 

rot {(pir)E) = ip{r) rot E + (p'{r)r'^RE , 

dw{(^{r)H)=(^{r)di\H + i^'{r-)r'^TH , (2.10) 
M((^(r)(E, H)) = ip{r)M{E, H) + ^'{r)r~^S{E, H) 

There is another correspondence between these operators. If we define the Fourier trans- 
formation 5" on g-forms in componentwise in Euclidean coordinates, then the mapping 
3^ : L^'"? — )■ L^''' is unitary and the well known formulas 

'J{dP E) = il"l Id" 5(E) , S'iE) = (- i)l°l3^(Id" E) 

and clearly 3^{AE) = —r'^J{E) , A3'(i?) = —3^{r'^E) hold for g-forms E . By elementary 
calculations we get 

g^rot = i/?3^ , ?div = ir3^ , ■JM = iS3^ , (2.11) 

rotJ=-i3^i? , divJ=-i3^T , M-J = -iS'S . (2.12) 



3 The time- harmonic problem 

Let r > and e , /i be two r-admissible transformations on q- resp. (g + l)-forms as 
well as M, A be as in (11. 7p . As mentioned above we want to treat the time-harmonic, 
inhomogeneous, anisotropic (generalized) Maxwell equation 

{M + iuK){E,H) = {F,G) 

with frequencies 

u eC+ := {z eC: \mz > 0} 
A substitution like x := ax , H := f3H allows us to suppose w. 1. o. g. 

£^0 = A^o = 1 and thus A = Id-|-A . (3.1) 

To shorten and simplify the formulas we always want to assume (13. ip throughout this 
paper. 

Now let us introduce our time-harmonic solution concept. From the skewadjointness 
of the two operators 

rot : R^iQ) C L2'''(fi) ^ L'^''^+\n) 
div : C L2''?+1(1]) ^ L^^'^in) 

to each other we obtain the selfadjointness of 
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with 

M{E, H):=i A-^M{E, H) = i{e-^ div H, ^"^ rot E) 
Here ,yL^'^(f2) := L^'''(f2) is equipped with the scalar product (z/ - , ■ )L2,q{n) ■ This suggests 

Definition 3.1 Let cu e C\R and (F, G) e L'^''^{n) x L^'i+\n) . Then {E, H) solves the 
problem Max(A, w, F, G) , if and only if 

(i) (F,i7) G R«(^]) X D''+l(^]) 

(ii) {M + iuA){E,H) = {F,G) 

The selfadjointness of M yields the unique solvability of Max(A, F, G) for each fre- 
quency ui e C\]R and all [F, G) G L^'''(f2) x L^'^+^(fi) . We denote the continuous solution 
operator by 

It can be seen easily that the spectrum of M is the entire real axis. Thus we expect, e.g. 
from Helmholtz' equation that we have to work in weighted L^-spaces and with radiating 
solutions to get a solution theory for real frequencies. 

Definition 3.2 Letco eR\ {0} and {F, G) G Lf^^(^]) x Lf^^+^(fi) . Then {E, H) solves 
the problem Max(A, cu, F, G) , if and only if 

(i) (F,ff)GR^^ ,(fi)xDf_\(fi) , 

^2 ^2 

(ii) {M + iuK){E,H) = {F,G) 

(iii) {r-^S + \d){E,H) G L'•^(^^) X L''«_+^(fi) 

> 2 ^2 

Remark 3.3 We call condition (iii) the 'Maxwell radiation condition' or 'radiation con- 
dition'. This condition generalizes the classical (N = 3 , q = 1) Silver-Mller incoming 
radiation condition for Maxwell equations (see (11.41) ) 

H - E eLl^^in) , ^x E + H eLl_^(n) 

We note that the radiation condition reads 

[r'^TH + E, r-^RE + H) e L^'*^ i (fi) x L^-'^Y {^) 

> 2 -^2 

Furthermore, we need 
Definition 3.4 We define 

P := {w G C \ {0} : Max(A, w, 0, 0) has a nontrivial solution.} 

and for u E C \ {0} 

l^{Max, A, u) := {{E,H) : {E,H) is a solution of Max(A, w, 0, 0) .} . 
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Clearly we have P C M\{0} and K(Max, A, = N{M-uj) = {(0, 0)} for w G C \ M . 

Similar arguments like those leading to the main result of the first part of |22] prove 
the following theorem. Therefore these do not have to be repeated here. We note that 
essentially we need two a priori estimates. Then the time-harmonic solutions are obtained 
by the limiting absorption principle. For details we refer the interested reader to fi3\ 
Kapitel 4] . 

Theorem 3.5 Let t > 1 andu eR\ {0} . 
(i) For allteW 

?^(Max,A,a;) = N{M-iu) 

c (R?(fi) n e-\Bt{n)) X n fi~\Rr\n)) , 

i.e. eigensolutions decay polynomially. 
Additionally let Q have the MLCP. Then: 
(ii) !)\f(Max, A, w) is finite dimensional. 

(ill) P has no accumulation point inM. \ {0} . 

(iv) For every {F,G) G L^'|(i7) x L^'|^^(i7) there exists a solution {E,H) of the problem 

^2 > 2 

Max(A, Lo, F, G) , if and only if 

A ((^'G'),(e,/^)>,,,(^)^,,,,,(^) = . (3.2) 

(e,/i)G>f(Max,A,a;) 

The solution can he chosen, such that 

(A(i5;,ff),(e,/^)>,,,(^)^,,,,,(^)=0 (3.3) 
holds for all (e, h) G K(Max, A, u) . By this condition {E, H) is uniquely determined. 

(v) The solution operator introduced in (iv), which we will denote by C^^ as well, maps 

(L2'9(fi) xL2''2+i(fi)) nX(Max,A,u;)^ to {Rt{Q) x nJ^{Max, A,u)^^ con- 

tinuously for all s, —t > 1/2 . 

Here we denote the orthogonality corresponding to the (A ■ , ■ )L2,g(j^)xL2.9+i(n)-scalar 
product by ±a and we put ± := ±id . Moreover, using the same technique introduced by 
Eidus in |1] for the classical Maxwell equations we get 

Corollary 3.6 Let t > 1, lu E R \ {0} and {E,H) G K(Max, A, w) . // additionally 
{e,fi) G C^'''(S) X C^''^^^(S) with bounded derivatives for some exterior domain E G Q , 
then 

exp(tr) ■ {E,H) G {R'^{Q) H e-^B%Q)) x {B''+\Q) n fi-^R''+\Q)) , 
exp(tr) ■ {E,H) G H^'^iE) x H^'^i+^E) 

hold for all t eM. and for all exterior domains S C 2 with dist(S, dE) > , i.e. eigenso- 
lutions decay exponentially. 
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Remark 3.7 The polynomial resp. exponential decay of eigensolutions holds for arbi- 
trary exterior domains , i.e. Vt does not need to have the MLCP. 

Remark 3.8 // the media are homogeneous and isotropic in the outside of some hall, 
i.e. supp A U (M^ \ VL) dUp for some p > , then 

supp(E,if) c nn Up 

for all u G M\ {0} and {E, H) G !N(Max, A, uj) , since in this case {E, H) solves Helmholtz' 
equation 

{A + u^){E,H) = {0,0) 

in Ap and therefore by Rellich's estimate {123] or [9, p. 59]) must vanish in Ap . If the 
principle of unique continuation holds for our Maxwell system, then 

K(Max,A,u;) = {(0,0)} 

Moreover, using the a priori estimate of the hmiting absorption principle and some 
indirect arguments followed by the (trivial) decomposition of Ll''^{Q) from [151 Lemma 
5.1] we are able to prove stronger estimates for the solution operator as the ones given 
in Theorem 13.51 (v) . 

Corollary 3.9 Let r > 1 , s, -t > 1/2 and K mC+\ {0} with Z n P = as well as 
have the MLCP. Then 

(i) there exist constants c > and t > —1/2, such that the estimate 

II ^'^(^'^)||R?(f7)xD?+i(Q) + II ^'S' + Id)£^(F,G')||L2,,(^)^L2'''+i(Q) 

holds true for all u e K and {F,G) G L^'^(i7) x L^'^+^(i7) . Especially the operator 

: Ll'%Q) X Ll'''+\Q) R?(r]) x D^'(^]) 
is equicontinuous w. r. t. uj a K ; 
(ii) the mapping 

C : K ^ fi(L2'9(l]) X L2'9+i(l]),R?(l]) X Df ^(fi)) 
UJ I — > 

is (uniformly) continuous. {Here we denote the bounded linear operators from some 
normed space X to some normed space Y by B{X, Y) .) 
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4 The static problem 

o 

To introduce our static solution concept we remind of the special forms 3*^(17) , 'B'^~^^(Q,) 
from [HI section 4] and the 'static Maxwell property' (SMP), which guarantees their 
existence and also implies the MLCP. (If Q is Lipschitz homeomorphic to a smooth exterior 
domain, then Q possesses the SMP.) To work with these forms we may assume that Q 
has the SMP, and restrict our considerations to ranks 1 < q < N . 

Definition 4.1 {E,H) is a solution of Max{A,0, f, F,G, gX,0 '^'^^^ data 

{f,F,G,g) G Lt-\n) X Ltlin) x if^f ^(fi) x Lf„f ^(^]) 
and {CO ^ ^ £_di+'^ ^ ^^ly 

{E,H) G (L^/^(fi)^RL(^])^5-lDL(^])) 

2 

^ 2 

solves the electro-magneto static system 

iotE = G , diYeE = f , 6^)L2,.(n) = , i = l,...,d'^ 

diYH = F , lot fiH = g , (/i//,6r')L2,,+i(n) = 6 , k = 1, . . . , d'^^' . 

Now we want to use [El Theorem 4.6] in the special case s = to solve the static 
problem Max(A, 0, /, F, G,gX,0 ■ t^iis let e , be r-C^-admissible with r > as well 

as 

where the latter is defined for g 7^ 1 , and for g 7^ 

Furthermore, for s = we put as usual o^'^i^) ■= oKi^) , o^'?(^) := o^;^(^) and 

W'?(fi) := wi{n) . 

Theorem 4.2 For every data (/, G, C) e W'(^]) and {F,g,0 e W«+i(n) i/iere exzsfo a 
unique solution 

{E,H) G (Rli(f]) n£-iDli(r])) X (Dlt'(f^) n/i-^fet'(^^)) 

of the electro-magneto static problem Max{A, 0, f, F,G, gX,0 ^'^^ corresponding so- 
lution operator is continuous. 
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Remark 4.3 For special data (0, G, 0) G W«(^]) , (F, 0, 0) G W+\n) , i.e. 

we will denote the corresponding continuous solution operator by 

Co : o3'{Q) X oK^+^(f^) ^ (fei(fi) X Dlt^(fi)) n A-i(o©li(l]) x ofct^(^)) ■ 
VKe note that Co even maps oD^(fi) x ot+H^) io 

(RUin) X D^+}(^])) n A-HoDLi(^^) x ofe(n)) 

continuously for all 1 — N/2 < s < N/2 . 

5 Low frequency asymptotics 

To approach the low frequency asymptotics of C^ we first have to be sure that P does 
not accumulate at zero. To this end first of all we derive a representation formula for the 
solutions of the homogeneous, isotropic whole space problem, i.e. Q := and A := Id, 
with the help of the fundamental solution of the scalar Helmholtz operator in 

A + u^ , wgC+\{0} 

This one can be written as 

^Lo,u{x) = ip:^^u{\x\) with ip^^i,(t) = CNOj''t~''Hl{ut) , 

where the constant cat only depends on the dimension N and Hl{z) represents Hankel's 
function of first kind for the index u := [N — 2)/2 . From now on we may additionally 
assume N to be odd, since then by the properties of Hankel's function (see e.g. [10] or 
[9|, p. 76]) v^^ y and its first derivative can be estimated by 

<c(t2-^ + t^) , <c(t^-^ + t^) (5.1) 

uniformly in t G IR+ and u & K ^ C+ with some constant c > depending only on N 
and K . 

From Remark 13.81 we have (in the case Vt = M^) 

?^(Max,Id,u;) = {(0,0)} 
Thus is well defined on the whole of L^'^ x L^'^^^ , if we denote C^ in the special case 

> 2 ^2 

= and A = Id by . Let w G C+ \ {0} and (F, G) G 0°°'^ x C°°''?+i . Looking at 
iE,H) ■.= L^{F,G) we get 

{E,H) G (H^'% n 0°°'") X (H^'^Y n C~'''+^) 
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by regularity, e.g. [HI Satz 3.6]. Applying (M - ico) to (M + iuj){E,H) = {F,G) and 
using iuj{dw E,Tot H) = (div F, rot G) we observe, that {E,H) satisfies 

{A + u^){E, H) = [M - iu - -n) {F, G) =: (/, g) G C°°'^ x C°°'^+^ (5.2) 



with □ := A - M2 



rot div 
div rot 



. We obtain {E, H) = (e, h) , where (e, h) is the 



unique radiating solution of the whole space problem 

{A + uj'){e,h) = {f,g) , 

(e, h) e Hj^ , X Hjn' , 
exp(- ia;r) ■ (e, h) e H^'S x H^n' 

^2 ^2 

For nonreal frequencies u G C+ \ M this is trivial, because again [131 Satz 3.6] yields 
{E, H) e H2'<? X H^-'J+i . But then (E, i/) = (e, /i) holds for real frequencies w G M \ {0} 
as well, since one receives the solutions of both radiating problems with the principle of 
limiting absorption. 

Using the representation formula for the solutions of the scalar Helmholtz equation, 
which, for instance, can be found in P pp. 78/79, Remark 4.28], we can represent the 
Euclidean components of our forms E = Ej dx^ and H = Hj dx'^ by 

Ei = fi-k , Hj = gj-k 

Here we denote the scalar convolution in by -k . For suitable g-forms e = cj dx^ and 
h = hi dx^ (Euclidean coordinates) we define the convolution 

e-k h := ei -k hj 

simply as the sum of the componentwise scalar convolutions. Furthermore, we have for 
suitable forms the rule of partial integration 

lot e -k h{x) = e -k diY h{x) . (5.3) 

With the special forms 

(J)-^ ■= $ . dx^ 

we get the representations 

Ei = f-^K,u , Hj = g-k^i^^ 

i.e. reminding of (15.20 

= (divG-iwF- -rotdivF) , (5.4) 

i/j = (rot F-iwG- -div rote) . (5.5) 



00 
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Our next goal is to use the partial integration formula f l5.3p to remove the second deriva- 
tives from F and G . Let us look at 

(divG)^<, 

for example. Because of the compact support of {F, G) we do not have to pay attention 
to the integrability of ^^^i, at infinity. By (15. ip we can estimate and V$aj,i/ in Ui by 
\^Lo,u\ < c ■ r^-^ , |V$^,^| < c ■ r^-^ and thus we have ^^^^ , V$^,j. G ^\Ui) . With the 
cut-off functions 

i^niy) ■= 'n{n ■ \x - y\) , neN , 
which satisfy |VV'n(2/)| < c ■ \x — y\^^ uniformly in n, we have 

Therefore (15. 3p yields 

(div Gn) = G„, -k rot 

with G„ := ipn ■ G and we obtain 

(divG)^$^,, = G^rot$:,, 

by passing to the limit n — ?■ oo and using Lebesgue's' dominated convergence theorem. 
Using these partial integrations in (15. 4p and (15.51) we finally get the representations 

Ej = Gk (rot $i J -iuF^ <l>i^^ - -^(div F) * (div J , (5.6) 

to 

Hj = F^ (div ^i^,) -icoG^ ^i^, - -{rot G) k (rot <^i^,) (5.7) 

UJ 

for any (F, G) G C~'« x C°°'«+i and {E, H) = L^{F, G) . 

Theorem 5.1 Let ^ co e K (£ C+ and s G (l/2,iV/2) , t := s - {N + l)/2 as well as 
{F, G) G X R^^^ . Then for {E, H) := L^{F, G) the representation formulas 

E={Gk (rot <|.^, J - i u;F ^ - - (div F) ^ (div $[,,)) ■ dx' , 
= (F * (div ^i^,) -icoGk ^i^^ - -{rot G) k (rot $i J) ■ dx' 

UJ 

hold in the sense of Lf''^ resp. Lf''^'^'^ . Furthermore, there exists a constant c > , such 
that 

+ ^ ■ ll(div^,rotG')||L2.,-i^L?-+0 
for allO^uj eK and (F, G) G x Rf+i . 



20 



Dirk Pauly 



Proof We choose a sequence C C°°''^ x 0°°'''+-^ converging to {F,G) in 

X R^^"*^ as n — 7- oo . Theorem 13.51 (v) yields the convergence of {En, Hn) := L^{Fn, Gn) 
to {E, H) G X D^+^ in R^^ x D^+^ since t < -1/2 . By ([51]) and ([EZD we may represent 
the forms {En,Hn) and observe that the involved convolution kernels essentially consist 
of Vu},u ° r and ^p'^^ o r . Using (15. ip these functions can be estimated by 

k-,. W I > W I < c (r'""^ + r^-^ + r^) < c (r^-^ + r^) 

uniformly in r and a; . From McOwen [H, Lemma 1] we obtain, that integral operators 
with kernels like \x — map continuously to Lf , if 

-N/2 <t<s<N/2 

As a direct consequence the right hand sides of (15.61) and (I5.7p define continuous linear 
operators from to . This proves the asserted representation formulas. By the differ- 
ential equation it is sufficient to estimate \\L^(F, G) |L 2,? . 2,9+1 • The uniform boundedness 
of the convolution operators w. r. t. 7^ a; G -ft' in the representation formulas yields the 
desired estimate, which completes the proof. □ 

For 7 G M+ we put 

C+_^ := {to G C+ : l^l < 7} 

From now on we assume that Q possesses the SMP, q ^ and e, /i are r-C^-admissible 
with order of decay 

r>(iV + l)/2 

We note that here it would be sufficient to demand the asymptotics 

i, fi,dni,dn fi = 0{r~'^) as r — )■ 00 , n = l,...,N 
Lemma 5.2 Let s G (1/2, N/2) and t := s - {N + l)/2 . 

(i) P does not accumulate at zero. In particular P has no accumulation point and there 
exists some u > , such that P fl C+^i, = . 

(ii) is well defined on the whole ci/L^'|(i7) x L^'|^^(r2) for all 00 G C^^^ \ {0} . 

> 2 ^2 

(ill) There exist constants c > and < u < u , such that the estimate 

II ^UliF,G)\\^2,g^^^^^2,g + l^^^ 

< C - (||(^,G')||L2.?(f^)xL?'«+^(n) + I'^l"^ ■ l|(^i^^'^°*'^)|lL?''-i(f7)xL?-'+2{C) 

+ ■ ^ |(F,6,«)L2„(n)| + ■ |(G,6f ^)L2„+i(n)|) 
1=1 i=\ 

holds true for all cu G C+,<i \ {0} and (F, G) G D^((]) x Rf+l(^^) . 
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(iv) Especially for all (F, G) G oD^(f^) x oK^+^(r]) and u E C+,^ \ {0} 

II >Co.(i^,G')||j^2,,(^)^L2-9+i(f7) ^ C - ||(^,G')||L2,,(f^)^L2''+i(n) 

The II ■ ||L2''j(Q)xL^'«+^(f7)"'^'^'"'^ ^^/'^ /ianc? sz(ies 0/ (iii) and (iv) may 6e replaced by 

the natural norm m [Til{VL) n e-^T>l{VL)) x (/i-iR^+^(fi) n Dt^+^(r])) . 

Proof First we prove tlie following: 

For all cl; > , s e (1/2, N/2) and t := s - {N + l)/2 there exist constants c, ^ > , 
such that the estimate 

IK^'^)|lL?''(n)xL?'«+i{n) 

< c- (||(^,G')||L2,,(f^^^L2>'+i(^) + ||(^,^)||L2,,(nn{/,)xL2.9+i(f7nc/,) (5.8) 
+ \uj\-' ■ ||(divF,rotG)||L2.-i(^^^)^L?.+2(^^^J 

holds for all oj G C+^^j \ {0} , all 

(F,G) G (L^'''(^]) n D^(A.J) X n R^'IAJ) 

and all solutions {E, H) of Max(A, uj, F, G) . 

Let {E, H) be a solution of Max(A, cj, F, G) and (F, if) the extension by zero of //(-E, H) 
to . This one satisfies the radiation condition, is an element of R'^ 1 x D"^ 1 , even 



of X H^^Y by [la Satz 3.6], and solves 

< 2 < 2 

{M + iu){E,H) = r]{F,G) + CM,r,{E,H)-iuA{E,H) =: {F,G) G x Rf ^ 

in since r > (A^ + l)/2 > s + 1/2 . Thus we^obtain {E, H) = L^{F, G) and Theorem 
15. II yields a constant c > independent of u , (F, G) or (F, if) with 

(F, ii) 2,q . 2,9+1 



< c ■ ( II G) II + ' ■ II (div F , rot G) 11^2.,-! ^^^'^+2 



(5.9) 



Furthermore, by the differential equations we get 

iwdiveF = divF , i a; rot /iii = rot G (5.10) 

in Arf^ and 

iwdivF = divF , iwrotF = rotG' (5.11) 
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in M . Combining (I5.9p and f lS.lip we have 

I|(^'^)|Il?-''(Q)xL2-''+1(Q) 

< c - (||(^,G)||l2.,(^)^l^'«+1(^) + ||(^,^)||L2,,(nni/.2)xL2.^+i(nni/.2) (5.12) 
and using flS.lOp we estimate the last term on the right hand side by 



I (div E, rot H) I 



- 2,5+2 



< C- -f^)||L2,9(Qn(7r2)xL2>9+i(f2nC/r2) 

+ II (div E, rot H) ||l^..-(s.pp,)xL^-+^(supp,), 

< C- (^||(-^; -f^)||L2,9(QnC/r2)xL2.'J+i(QnC/r2) 

+ II (div rot fiH) ||l2,.-1(,,pp^)^l?-+^supp,) 
+ ■ ||(divF,rotG)||j^.„_.(^^^)^L.„+.(^^^)) 

< C- (^||(-£^; -f^)||L2,9(f2nC/^2)xL2.'J+i(f2nC/r2) 

+ II ^)|lH,^L^(supp^)xHiL'^'(suppr,) 

+ \u;\-' ■ ||(divF,rotG)||L..,_.(^^^^^L..,+.(^^^J . 

Inserting this estimate into fl5.12p . using the regularity result [I3l KoroUar 3.8 (i)], the 
differential equation as well as fIS.lOp we finally get 

l|(^'^)|lL?-9(f7)xL?-9+i(n) 

< c- (||(^,G)||l2.,(j^)^l2'9+1(q) + U^^^^Wh^'XmxL^'^+'in) 

+ ■ ||(divF,rotG)||L2„-.(^^^^^L2„+2(^^^J . 

By r > (iV + l)/2 we have s - r <t and thus flS.Sp follows. 

If we now assume that is an accumulation point of P or the estimate in (iii) is false, 
then there would exist a sequence {uJn)n£N C C+ \ {0} tending to zero and a data sequence 

((F„,GO)„,i, ^ (D^^) X k^\^)) nK(Max,A,a;0^ 
as well as a sequence of normed solutions {En,Hn) to (M + iijJnA){En, Hn) = {Fn,Gn) 
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with ||(^n,i^n)||L2.,(^)^L?-+^(n) = ^ 

l^nl"^ ■ ||(divi^n,rotG'„)||L2,g-i(f^)^L?'''+^(f7) 

|a;„|-i-|(G„,6f^)L2„+i(f,)| ^0 , k = 1, . . . , d'^+' 

(in the case of (iii) we have of course {En,Hn) = By the differential 

equation we get ici;„(dive_E'„, rot = (div rot G^) and thus 

+ ||(div£E„,rot/ii/„)||L2,..i(^^^L^.«+2(^) 

Consequently {En, Hn) is bounded in 

and thus the MLCP yields a subsequence, which we also denote by , con- 

verging for every t < t in L?''^(r2) x L?''^^^(r2) . Because of (15. 13^ this sequence even 
converges in 

(R?(f]) ^£-lD?(^])) X (/i-^Rf'(fi)nD?+'(n)) 

to the Dirichlet forms, let us say 

Since t = s — {N + l)/2 G (— A^/2, —1/2) we may assume w. 1. o. g. t > —N/2 . Therefore 
by (12.81) we obtain 

{E, H) e s^'i^) X fi-\-i:K''+\Q) 
For i = 1, . . . ,d'^ we compute 

\l^n\ ■ I V-Tn, 0£/L2.'?(f1) I > U 

o o 

= \i^n\~ ■ (div/7„,6^)L2,9(n) + iu;„(£:i?„,6^)L2,q(n) 

=0 

= |(£^-En,&£)L2,9(Q)| |(£^-£^,&f)L2>9(Q)| , 

o 

i.e. E G B^(r2)-^^ . Analogously we see H e B'?+^(n)^''. Thus (S,i/) must vanish and 
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finally fl5.8p yields constants c, g > independent of n with 

1 = WiEn, Hr,)\\^2,,^^^^^2,,+l 



<c - 



l2'9(Q)xL2'9+1(Q) 



+ \un\ ^ ■ ||(divF„,rotG„)||L2.,-i(^^^^)^L2'''+2(A.o) 
+ IK-^"' -^")|lL2.9(nnc/e)xL2.9+i(Qnc/p)) ' 
a contradiction. □ 

We are ready to prove our main result: 

Theorem 5.3 Let s G {1/2, N/2), t := s - (A^ + l)/2 and Co he from Lemma [Qi 
Furthermore, let {uJn)nim C ^+,Cj \ {0} 6e a sequence tending to and 

6e a dato sequence, such that 

-iu;-i(divF„,rotG„)^ (/,(?) m 

-iuj~'^{Fn,hl)i2,q(^^^ ^^^—^ Qi zn C , £=l,...,d'^ , 

-ia;.^^(G'„,6^^^)L2.<j+i(f7) ^^-^ ^fc m C , /c = 1, . . . , ci^+^ 

/ioW. Then {{E„,, H„))^^^^ := ( £^„(F„, G'„))^^pj converges for all i < t in 

to {E, H) , the unique solution of the static problem Max(A, 0, /, F, G,g,(,0 ■ 
Proof From Lemma [5.21 we get the boundedness of in 

Thus by the MLCP we can extract a subsequence, which we will denote by [{En, Hn))_^^-^ 
as well, such that 

Hn) {E, H) in h^in) x Lp+^(n) 

holds for all t G {—N/2,t) . The differential equation 
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and the assumptions yield 

M(E„, H^) ^ {F, G) in X L?'^+^(fi) , 

(div£E„,rot/i//„)^^ in . 

For /c = 1, . . . , (P^^ we compute 

6|"^^)l2,9+1(q) = (roti?„,6^"^^)L2,g+i(n) 6^"^^)L2.<7+i(n) "~^°°> ^fc 

l^n ^ -V ' '^ra 

=0 

o n— s^oo ~ ~ 

and analogously {sEn, > Q ioi I = 1, . . . ,df . Thus (-E, if) is an element of 

^2 ^2 ^2 ^2 

solving the electro-magneto static system 

rotE = G , div^ = F 

div eE = f , rot fiH = g , 

For the difference (e, h) := {E, H) — {E, H) we obtain 

{e,h) G {,:Kl ,^{Q)nB''{Q)^^) X {fi-\-i:Kl^\{Q) nB''+\Q)^-) 

^2 ^2 

and even (e, /;.) G L^''(f2) x L^''^'*"^(r2) again by f l2.8p . Thus (e, /i) must vanish and because 
of the uniqueness of the limit {E,H) = {E,H) even the whole sequence iin))„gpj 
must converge to {E, H) in L^^(^]) x L^^+^(^]) . □ 

Corollary 5.4 Let s , t , Co he as in Theorem{5lE and {F,G) G o'^ii^) x ofe+H^) ■ 
r/ien the solutions Ci^{F,G) of the time-harmonic problem Max{A, u, F,G) converge for 

all t < t in R?(f^) x D?"^ (fi) to Cq{F,G) , the unique solution of the static problem 
Max(A, 0, 0, F, G, 0, 0, 0) , as cu e C+,i \ {0} tends to zero. 

By a similar indirect argument (see [131 Korollar 7.5]) we obtain 

Corollary 5.5 Let s G (1/2, N/2) , t := s - {N + l)/2 , u be from Lemma\5^ and Bs,t 
be the Banach space of bounded linear operators from the Hilbert spaces 

oBl{Q) X oK^+^(^^) to R1{Q) X 

Then \\ £^ Ws^-t is uniformly bounded w. r. t. u E C+,(2, (even for u = !). Moreover, the 
mapping 

C : C+,<2, — > B^i 

oj I — > C,,, 



is (uniformly) continuous for all t < t . 
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Remark 5.6 Clearly Rt{n) x Bf^^Q) resp. R?(f2) X D?'*'^(f2) may be replaced by its 
closed subspace 

(R?(n) X n K-\om^) X 

resp. 

{iq{Q) X Df ^(fi)) nA-^(oD?(^^) X 

Corollary 5.7 Let s , t , uj , {ujn) be as in Theorem \5.3\ as well as 

which may be decomposed by /75|. Theorem 3.2 (iv)], such that 
{Fn, Gn) = A(F^, G^) + (F^, G^) 

with 

{F:,G'J e (oM^(^^) +LmB'^(fi)) x (oDf ^(fi) + LinB^+i(fi)) , 

Moreover, let ((F^jGiJ^)) converge to some {Fd,Gr) in L^'''(fi) x L^'^+^(r2) as well as 
( - cont;er(/e to some (E^ if'^) m L^'^(fi) x L^'^+^(l]) /or all i < t . Th 



en 



orm 



{{En,Hn)) := {C^„{Fn,Gn)) convcrgcs for all i < t m L^-'^(fi) x L^-'^+X^) ^ the f( 
{E,H) = {E^,H'') + Co{F,,Gr). 

Proof converges to Co{Fd,Gr) by Corollary 15.41 Moreover, of course 

holds. □ 



6 Inhomogeneous boundary data 

We want to finish this paper by discussing inhomogeneous boundary data. 

Recently Week showed in [28j, how one may obtain traces of differential forms on 
Lipschitz boundaries. Let H be a bounded Lipschitz domain in . Then we know from 
[28| Theorem 3] the existence of a linear and continuous tangential trace operator (using 
for a moment the notations from there) 



r : R^(S) 
E 



L*E 
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Moreover, he proved in Theorem 4] that T is surjective, i.e. the existence of a 
corresponding hnear and continuous tangential extension operator (a right inverse) 

Let e be a 0-admissible transformation. Applying the usual Helmholtz decomposition 

L2'''(S) = rot R9-i(S) ®e £^"(2) ®e e-^diY D5+i(S) 
we receive a linear and continuous tangential extension operator with range in 

n e-Miv D9+i(2) C R^(S) n e-\B''{E) 

If we assume now that Q possesses a Lipschitz boundary (This implies the SMP.), then 
we get by an usual cut-off-technique for any s G M a linear and continuous tangential 
trace operator 

7, : RUn) — > widn) := R-^/^'''{dn) 

and a corresponding linear and continuous tangential extension operator 

o 

satisfying 7^-7,- = Id on TZ^{dVL) . We note that the kernel of 7t- equals Rs(f2) and that 7,- 
may be defined even on R^Q^(f2) . 7^ may be chosen, such that supply,- A C fl Ur^ holds 
for all \en\dVt). 

Let (F,G) G Lf^^(n) X and A G TV{dQ) be some boundary data. We want 

to discuss the solvability of the time-harmonic Maxwell system 

{M + iujK){E,H) = {F,G) , irE = \ (6.1) 

using the results obtained so far. By definition we have 

EA:=7.AGRU^)n5-iDUfi) 

and with the ansatz 

(E,iJ):=(E,^) + (EA,0) (6.2) 

the equations fl6.ll) turn to 

(M + iwA)(^,if) = (F,G) , 7.^ = (6.3) 

_ _ _ o 

with (-F, G) := {F,G) — {i us Ex, rot Ex) ■ Thus we are looking for E G RfQ^(fi) and we 
can use the results from the previous sections. Moreover, for any s G M we clearly have 

(F, G) G L^'^n) X Ll'^^\Q) ^ (F, G) G L^'^Q) x L^'^^^Q) 



28 



Dirk Pauly 



and for nonreal frequencies u; G C \ R and {F,G) G L^''^(f2) x L^''^+^(fi) we easily get 
unique square integrable time-harmonic solutions 

{E,H) :=£^(F,G) + (Ea,0) 

= C^{F,G) - C^{iue%X, rot %X) + {%X,0) G R'^iQ) x B'^+\Q) . 

We denote the continuous solution operator by 

: L^''i{Q) xL'^'i+^Q) xTZ^idQ) — > R'^ {Q) x W^+^Q) 
(F,G,A) ^ iE,H) 

and note £^ = iSt^( ■ , ■ , 0) . 

To establish a solution theory for non vanishing real frequencies a; G M \ {0} and data 
{F,G) G L^1(n) X L^1+^(n) with our Fredholm theory from Theorem 13.51 we consider 

> 2 > 2 

r-admissible transformations {e, /i) with some r > 1 . Using the ansatz (16. 2p we only have 
to guarantee 

± X(Max,A,tj) 
Let (e, /i) G 3\f(Max, A, u) . We compute 

= ((F,G),(e,/i)>L,, 

9(n)xL2.'J+i(n) ~ (rot-E'A, ^)l2.9+1(Q) + (i?A, il^^e)L2,g(n) 

= ((^^G*), (e,/i))L2,,(j^)xL2,.+i(n) ~ (r0tEA,/i)L2..+i(f7) - (^A,div/l)L2,g(n) 

= ((-^><^), (e, /i))L2,9(j^)xL2.9+i(n) ~ (^t-^A,^n/i)L2,^r(Q)xL2.<7+i(f1) 

with Tt^ := ($, rot $) and T„^ := (div ^, ^) . 



Remark 6.1 Assuming more regularity of Q , i.e. Q E , and fi, i.e. ^ E , hy 
Stokes' theorem 

(TtEA,T„/i)L2,<j(n)xL2.^+i(n) = (7r^A,7n^)H-i'j(an) " ^"^^ H^i^ca n) 

holds, since then by regularity h is an element of and thus 'jnh is an element 

o/ H^/^'^(9fi) , where 7„ = ± ® l** denotes the usual normal trace. {^Here ® denotes 
the star-operator on the submanifold dQ of Q and (■, ' )jj-^,9(-gQ-) duality between 

U-h'i{dn) andU^'''{dn) .) 

These considerations yield the following solution concept for w G M \ {0} : 

We call {E, H) a solution of the radiation problem Max(A, u, F, G, A) , if and only if 

. {E,H)eRl_,{n)x-Dl+\{Q) , 

^2 ^2 
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• {M + iujA){E,H) = {F,G) and 7,^ = A , 

• {r-'S + ld){E,H) e L^'\{n) X L^^'^\\n) . 

> 2 > 2 

Theorem 6.2 Lei (£:,/i) be r-admissible with r > 1 . For all tu E M. \ {0} , A G TZ'^^dVl) 

and G Lj|(fi) x LM^^(fi) t/iere ea;zsfo a solution {E,H) 0/ Max(A, F, G, A) , z/ 

-^2 > 2 

and on/y z/ 

/or a// (e, /;.) G !N(Max, A, cu) . T/ie solution can be chosen, such that 

{E,H) ±A K(Max,A,a;) 
Then by this condition the solution (F, H) is uniquely determined and the solution operator 

: L^'lin) X LH'^^n) X w'idn) x d'^+_'i(^) 

^2 ^2 ^2 ^2 , 

(F,G,A) ^ (F,if) 

where {E,H) = C^{F,G) — £(^(ia;e7^A,rot7^A) + (7^A,0) , is continuous in the sense of 
Theorem \3.5\ (v). 

Now we need an adequate static solution theory to describe the asymptotic behaviour 
of . 

We call (F, H) a solution of Max(A, 0, /, F, G, g, (, ^, A, x) , if and only if 
{E,H) G (R'^ ^{Q)ne-'Dl ^{Q)) x (^-^Rf ^^(^]) n Df ^^(^])) 

-^2 ^2 ^2 -^2 

and 

rot F = G , div F" = F 

div eE = f , rot fiH = g , 

7rF = A , 'jrfJ'H = >c 

hold. 

For the rest of this paper let g 7^ . From [HI Theorem 6.1, Remark 6.2] (in the 
special case s = 0) we get 

Theorem 6.3 Let {e,fi) be t-C^- admissible with r > 0. Then for all f G qID''^ "'^(^) ; 
F G oD^(fi) , C e ^ G C^'^' an(^ all G G oR'+H^) . 9 ^ oR^+^(^^) , A G 7^'?(9^]) , 
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X G TZ'^^^{dri) satisfying 

RotA = 7^G A yA^ (C, 6)L2,g+i(n) = (^t7TA,T„6)L2,q(Q)xL2.?+l(f2) 

= (rot7^A,6)L2,g+i(f7) 

Rotx = 7^5( A f\ (fi',6)L2.'j+2(f7) = (^t7r>^,7'„6)L2,<7+i(n)xL2.9+2(n) 

fceB9+2(c) 

= (rot7^x, 6)L2,g+2(n) 

t/iere exists a unique solution 

{E,H) e (Rli(fi) ^£-^Dll(^])) X (^-iRiV(^) n D!+^(^])) 

0/ Max(A, 0, /, F, G, g, C,, A, x) . The solution depends continuously on the data. 
Remark 6.4 Once again assuming more regularity of Q , i.e. G , we have 

(Tf7^A,T„6)L2,9(n)xL2.<7+i(Q) = ('^'7n&)jj-i,g^g^^ 

resp. 

(Tt7^x, r„6)L2,5+i(Q)xL2.'?+2{n) = (^'T'"^)H-2'9+i(an) 

Finally we are ready to prove our last result: 

Theorem 6.5 Let{e,n) be r-C^- admissible with t > {N + l)/2. Let s e {1/2, N/2) and 
t := s — {N + l)/2 as well as u be from Lemma lSTR Moreover, let {um)mem C C+_(i \ {0} 
be a sequence tending to zero and 

be some data sequences with 

IrGm = Rot Am , 

such that 

{Fm,Gm)^^{F,G) m L^.«(l]) X L2'«+1(1]) 
-iu;-i(divF„,rotG„)^^(/,^7) in ^ 

o m— >oo 

-'^^^{Fm^Kli^^'^ip.) in C , i=l,...,d'^ 



-(rot7^Am,6^"^^)L2.g+i(n)) '"^°°> 6 in C , k = 1, . . . ,d'^^^ 
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hold. Then {{Em, Hm))^^^ := (>S^„(F„, G^, -^m))^^^ converges for allt <t in 

(R?(i]) ^e-lD?(^])) X (/i-lRf^(^]) n Bf\n)) 

to {E, H) , the unique solution of the static problem Max(A, 0, /, F, G, g, (, ^, A, 0) . 
Proof From Theorem lO and flOj) we have {Em,Hm) = {Em,Hm) + (^a„,0) with 

^Xm ■= Ir^m, {Em,Hm) ■= Cui^{Fm,Gm) and 

Fm '■= Fm — iUJmSEx^ , Gm := Gm — ^Ot i^A^ 

Because of the compact support of Ex^ and the continuity of % we have 

Ex^ Ex := irX in R^(l]) n e-'Dl{n) 

for all s G M . Moreover, (F^, G^) fulfills the assumptions of Theorem 15.31 Thus {Em, Hm) 
converges for alH < t in 

{Yq{n)r\e-^T>'i{n)) x {^i-^-kf\n) n D?+'(n)) 

to (F, H) , the unique solution of Max(A, 0, /, F, G, g,(,^) with F = F,g = g,$,=^ and 

G = G-TotEx , f = f-dweEx , C = C - [(^^^a, &^)L2.,(n)] ti • 

We obtain {Em, Hm) "^°^> {E, H) := {E, H) + {Ex, 0) with the asserted mode of conver- 
gence and clearly {E,H) is the unique solution of the static problem 

Max{A,0,f,F,G,g,C,^,X,0) , 

which completes the proof. □ 
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